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Abstract 

A first-order model for a stock market assigns to each stock a return parameter and a 
variance parameter that depend only on the rank of the stock. A second-order model assigns 
these parameters based on both the rank and the name of the stock. First- and second-order 
models exhibit stability properties that make them appropriate as a backdrop for the analysis 
of the idiosyncratic behavior of individual stocks. Methods for the estimation of the parameters 
of second-order models are developed in this paper. 
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1 Introduction 

First-order and second-order stock market models are relatively simple stochastic models that man- 
ifest some of the stability properties of actual stock market behavior. These models arc descriptive 
as opposed to normative, and are constructed using data analysis based on actual stock markets. 
First-order models are stock-market models where the parameters for return and volatility are based 
on the ranks of the stocks. These models were introduced in Fernholz (2002) and developed in Ban- 
ner, Fernholz, and Karatzas (2005), and reflect the actual rank-based growth rates and variances of 
the stocks in the market. First-order models are asymptotically stable, and accurately reproduce the 
long-term characteristics of the market's capital distribution. However, these models are ergodic in 
the sense that each stock asymptotically spends equal average time at each rank, and this ergodicity 
property does not seem to be present in actual markets. This lack of verisimilitude is the motivation 
to consider the next level of complexity: second-order models. 

Second-order models are a form of hybrid Atlas models, where the return and volatility param- 
eters are based on the rank and the name (or index) of the stocks (see Ichiba et al. (2011)). While 
these models retain many of the characteristics of first-order models, the above ergodicity property 

is no longer present, and this produces a more realistic representation of actual stock market behav- 
ior. In second-order models, larger stocks tend to remain asymptotically among larger stocks, and 
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smaller stocks tend to remain among smaller stocks. This behavior is closer to that of actual stock 
markets, so second-order models provide a more accurate descriptive representation of stock market 
behavior. 

Estimation of the parameters for first-order models is fairly straightforward, and can be accom- 
plished without great ado. Second-order parameter estimation is somewhat more complicated. Here 
we shall focus on the growth-rate parameters, and find it necessary to rely on implicit methods to 
determine values for these parameters. Our purpose here is to develop techniques for estimating 
second-order growth-rate parameters, not to carry out an exhaustive examination of these parame- 
ters for an entire stock market. First, let us establish some formal definitions. 

A m,a.rket is a family of stocks X = (X\, . . . , X„) whose capitalizations are modeled by continuous, 
positive semimartingales that satisfy 

d 

dlogXi{t) = Gi{t) dt + Y. ^-^-W' (1-1) 

for i G M, where n < d, B = {Bi, . . . , Bj) is an M"^- valued Brownian motion defined on M, and the 
Gi and Siv are progressively measurable with respect to the Brownian filtration, with Gi locally 
integrable and Si^ locally square-integrable. The reason we define these processes on M is that in 
practice we are confronted with time series over a given block of time, and the analysis of these 
series can be performed in both forward and reversed time. Hence, we see a sample in time of the 
processes Xi, . . . , and draw our conclusions from this sample. 

Wc shall assume that for any i ^ j, the intersection sets {t : Xi{t) — Xj{t)} have Lebesgue 
measure zero, almost surely, and we shall also assume that there are no triple points, i.e., ifi < j < k 
then there is almost surely no t e K such that Xi{t) = Xj{t) = Xk{t). The general setting for our 
model can be found in Fernholz (2002) and Fernholz and Karatzas (2009). 

The value Xi{t) of the stock X^ at time t represents the total capitalization of the company at 
that time. If we let Z represent the total capitalization of the market, then 

Z{t)^Xi{t) + ---+Xn{t), 

and we can define the market portfolio to be the portfolio ^ with weight processes given by the 
market weights 

Mi(*)=^^^, tor 1 = 1,..., n. 

We shall assume that the market weight process n = (/xi, . . . ,Hn) has a stable, or steady-state, 
distribution, and that the system is in that stable distribution. We shall be interested in the relative 
behavior of the log-capitalizations or log-weights. If iJ.{t) is in its steady-state distribution, then the 
log-difference processes defined by 

log Xi{t) - log Xj{t) = log fii{t) - log/ij(t), 

for i,j = l,...,n, will also be in their steady-state distribution. 

Consider the ranked capitalization processes corresponding to the Xi (t) in descending order 

> ••• >^(„)W, 

and the corresponding ranked market weights 

M(i)(i) > ••• > 
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Let rt{i) represent the rank of Xi{t), and let pt be the inverse permutation of Vt (with ties in rank 
settled by the order of the indices), so 

Xi{t) = and X(^h){t) = Xp^(^k)(t), 

and, similarly 

MiW = and IJ.{k){t) = l^pt{k){t)- 
Hence, pt(fc) represents the index, or name, of the stock occupying rank k at time t. 

The ranked market weights . . . = (Mpt(i)(0) • • • )Mpt(n)(0) comprise the capital 

distribution curve of the market at time t. The capital distribution curves over several decades of the 
20th century can be seen in Figure 1, a version of which appears in Fernholz (2002). The curves in 
Figure 1 show the ranked market weights on December 31 of the years 1929, 1939, 1949, 1959, 1969, 
1979, 1989, and 1999. During that period, the number of stocks in the market increased over each 
decade, so the decade associated with each curve is clear from the chart. We see that the capital 
distribution curve of the market shows a certain stability over time, so the assumption that fx is in 
its steady state distribution would seem to be consistent with the observed data. 




o 



1 5 10 50 100 500 1000 5000 

RANK 

Figure 1: Capital distribution of the U.S. market: 1929-1999. 
The curves show the ranked weights at the end of each decade. 
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2 First-order models 



A first-order model is a stocic-market model in which each stock has constant growth and variance 
parameters that depend only on the rank of the stock by market capitalization. These models were 
developed in Fernholz (2002) and Banner at al. (2005), and can be constructed to reflect certain 
properties of actual stock markets. A first-order model that is based on an actual market will have 
a steady-state capital distribution curve that is about the same as the capital distribution curve for 
the actual market (see Fernholz (2002), Figure 5.6). 

A first- order model is defined by a system X = (Xi, . . . , X„) of the form 

d\ogXi{t) = gf:^(^i) dt + (7f^(j) dWi{t), 

71 n 

= ^gklhii) = kdt + ^C7felft(i) = kdWi{t), 
fe=i fe=i 

for i = 1, . . . , n, where gi, ■ ■ ■ ,gn are real constants, cii, . . . , cr„ are positive constants, and {Wi, . . . , Wn) 
is an ]R"-valued Brownian motion, and where ft{i) represents the rank of Xi{t) (analogously to rt{i) 
for the rank of Xi{t)). We shall assume that the gk satisfy 

gi-\ h 5n = 0, 

and 

m 

fe=l 

for m <n. With these parameters, the Xi form an asymptotically stable system, which means that 
the market weights 'jliit) = Xi(t) / {^Xi(t) H -|- X„(f)) satisfy 

lim t~^ log 'fiiif) =0, for z = 1, . . . , n, 

t— >oo 

and the limits corresponding to (2.1) and (2.3) below exist (see also Fernholz (2002), Definition 5.3.1). 

Suppose we have a market X, and suppose that its market portfolio /i is in the steady-state 
distribution. We define the asymptotic rank-based relative variances for the market by 

lim t-i(logM(fc))(t), (2.1) 

and the asymptotic rank-based relative growth rates by 

gk= lim 1 / Vlrt(i) = fed log Mi (i), (2.2) 
Jo 

and suppose these limits exist almost surely. Since these parameters are based on the market weight 
processes /Xj, they represent values relative to the market portfolio /x. 

For k < £, let Ak,£ be the local time of the nonnegative semimartingale log(/^(fe)//X(^)) > at 
the origin, and set AO, 1 = = An, n + 1. Since we have assumed that the X^ almost surely have no 
triple points, it follows that for £ > k + 1, the local time by Ak,£ is identically zero, so here we need 
to consider only local times of the form AA:, fc + 1, and we have 

rflog/i(fc)(t) = ^lrt(i) = kdlogiJ.i{t) + -dAk,k-\- l{t) - -dAk- l,k{t),a.s. 
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For k = l,...,n — 1, we can define the asymptotic local time 

\k,k+i= lim t-^Afc, k + l{t), (2.3) 

which exists almost surely, and define Ao,i = = Xn,n+i- It turns out that the estimation of the 
\k,k+-L is not difficult, and the procedure is described in the appendix of Fernholz (2002). It can be 
shown (c.f. Proposition 5.3.2 in Fernholz (2002)) that 

9fe = ^ (-^fc-i,fe - Afe,fc+i) , a.s. (2.4) 
holds for fc = 1, . . . , n, and it follows that gi H + 9n = 0. 

The smoothed values of cr^ and Qk for the largest 5120 stocks in the U.S. market for the decade 
1990-1999 are shown in Figures 2 and 3. Since the number of stocks in the market changed during 

the decade of 1990 1999, wc limit our attention here to the largest 5120 stocks, which is fewer than 
the number of stocks in the market at any time during that decade. The vahics in Figure 3 do not 
add up to zero, since the largest 5120 stocks arc a strict subset of the larger market. 

The first-order model X = {Xi, . . . . Xn) such that 

dlogXi{t) = gf^(i)dt + crf^(^i)dWi{t), 

n n 

= "^gklhii) = kdt + ^<7fclft(i) = kdWi{t), 

k=l k=l 

where ft{i) is the rank oi Xi{t) at time t, is called the first- order model for the market X. As we have 
seen, the growth and variance parameters for X are derived from the relative growth and variance 
parameters corresponding to the market weight processes jUj, not directly from the capitalization 
processes X^. 
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Figure 2: Smoothed values of cr^. A: = 1, . . . , 5120, for U.S. market: 1990-1999. 
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Figure 3: Smoothed values oi gu, k = 1, . . . , 5120, for U.S. market: 1990-1999. 
The Qk satisfy J2^=i9k ~ ^' where n=7000. 



First-order models are ergodic in the sense that 



lim 1 / Ift(i) = kdt = lim 1 / lXi{t) = Xrk){t)dt = -,a.s. 



T-yoo 



T->-oo 



(2.5) 



This ergodicity property does not seem to be present in real markets, but instead there exist asymp- 
totic occupation rates defined by 



eki= lim 1 / lrt{i) = kdt = lim 1 / lXi{t) = X,k){t)dt,a.s. 



(2.6) 



Here 9ki represents the fraction of time that Xi spends in the kth rank. The n x n matrix 9 = {9ki) 
is bistochastic, and we shall assume that all the entries are positive. For a first-order model, (2.5) 
implies that 9ki = 1/n for all i and fc, and since this docs not seem to characterize the behavior of 
real markets, we shall now consider a more general class of models. 



3 Second-order models 



A second-order model is a stock-market model in which each stock has constant growth and variance 
parameters that depend on the rank and name, or index, of the stock. Second-order models are 
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examples of hybrid (Atlas) models, which were discussed in Ichiba et al. (2011). A second-order 
model is defined by a system X = {Xi, . . . , X„) of the form 

dlogXi{t) = (7i + gft{i))dt + ai^ft{i) dWi{t) (3.1) 

n 71 

= (li + ^9klrt{i) = kjdt + J^'^iklhii) = kdWi{t), 
fe=i fe=i 

for i = 1, . . . , n, with constants gk, 7j and aik > 0, for i, A; = 1, . . . , n, and a Brownian motion W. 
In order for the Xi to be asymptotically stable, these parameters must satisfy 

5i H H 5™ = = 7i H H 7n, 

and, for any permutation tt G S„, 

m 

^{9k + l^{k)) < 0, for m<n. 
fe=i 

Here we are interested in estimating the growth-rate parameters 7^ and gk- For simplicity, we shall 
consider only rank-based variances, and assume that u?^ = cr^ for all i and k. 

It was shown in Ichiba et al. (2011) that a second-order model of the form (3.1) is asymptotically 
stable, and the asymptotic occupation rates 



A 



9ki= lim 1 

T->oo 



/ lft(i) = kdt (3.2) 
Jo 



arc defined for all i and k, almost surely. The matrix 6 = (6ki), like 6 in (2.6), will be bistochastic 
with positive entries. We can generate the first-order parameters and for X as in (2.1) and 
(2.2), with 

CT^= lim t-^(log/I(fe))(i), 

and 

/■T " 

gk= lim 1 / V'lft(i) = fed log /ij (i) . 

Jo tt 

With these parameters, it was shown in Ichiba et al. (2011) that, almost surely, 

n 

9k = gk + ^Okai (3.3) 

n 

= 7i + ^^fei5fe- (3.4) 
fe=i 

In matrix form, this can be expressed 

9 = 5 + ^7 
= 7 + ^5, 

where 7, g, and g are column vectors. Prom this we see that 

7 = -^5, (3.5) 

so 

g={ln-90^)g. (3.6) 
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4 Estimation of second-order parameters 



The first-order growth parameters gk for the market X can be estimated directly from the stock 
return time series; however, second-order growth parameters will have to be estimated indirectly. 
We wish to construct a second-order model that has first-order growth parameters equal to those 
of the market, and an occupation-rate matrix equal to the occupation-rate matrix of the market. 
Under these circumstances, as in (3.6), we have 

g={ln- e6'^)g, (4.1) 

and we wish to solve this equation for g, the vector of name-based growth parameters for the second- 
order model of the market X. If we can solve (4.1) for this g, then we can use (3.5), in the form 

7 = -e^9, (4.2) 

to generate the name-based growth parameters 7i for this second-order model. Let us first consider 
the matrix 6. 

The matrix 9 is bistochastic and we have assumed that all its entries are positive, so this also 
holds for 6^ and 96^ . By the Perron-Frobenius theorem (see Perron (1907)), the symmetric matrix 
96^ will have a simple eigenvalue equal to 1 with eigenvector e\ = (1,1,..., 1)', and all the other 
eigenvalues will have absolute value less than 1. Hence, /„ — 99"'" has rank n — 1 and its kernel is 
generated by ei, so the condition that the gj- sum to zero means that g is orthogonal to this kernel, 
and this ensures a unique solution to (4.1). 

Unfortunately, it seems to be essentially impossible to estimate 9 with any reasonable accuracy, 
so although we can use this matrix to prove the existence and uniqueness of g, in practice we cannot 
actually solve equation (4.1). Instead, let us consider (3.3) in the form 

n 

gk=9k + ^Okili- (4.3) 

i=l 

We can use this equation to generate the g^ recursively, and then estimate the 7j from the returns 
data and the gk- 

Let us assume that the market X is defined for all f € M, that the weight process n for X 
has a stable distribution, and that is in that stable distribution. We can then define the time- 
reversed market X with stock capitalizations Xi{t)=Xi{—t) and weights 'jli{t)=fii{—t), and with 
this definition we can define the expected backward occupation rates similarly to (2.6). Since the 
weight process is in its steady-state distribution, the limits of (2.6) will be the same at plus and 
minus infinity, so the forward and backward expected occupation rates 9ki will be equal. The results 
of Bertoin (1987) imply that the forward and backward asymptotic local times Afc, A; -|- 1 will also bo 
the same, so the forward and backward versions of the are equal. Hence, it follows from (2.4) that 
the forward and backward gk are equal. In this case, (4.1) implies that the forward and backward 
values of the gk are equal, and from (4.2), we see that the forward and backward 7i are also equal. 
Quadratic variation is invariant under time reversal, so the forward and backward ak will be the 
same. Hence, the first- and second-order models for X are the same as the corresponding models 
for X, and this allows us to use both X and X to estimate the second-order parameters. 

In order to estimate the second-order parameters, it is necessary to observe the movement of 
market weights forward and backward in time. To this end, we define the concept of flow in a 
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market. The forward flow cjik of the market at rank k is defined for r > by 

^.(r)Alin.l/^og('^^fi^)di, 

T-).<x) Jo V M(fe)W ^ 

and the backward flow <pk of the market is defined by 

^.(r)-liml/^og(^^^4^)d. 

In Figure 4 we see the exponential of forward and backward flow for the largest 250 stocks in the 
U.S. market over the decade from 1990 to 1999. The plots show the average exponential flow of each 
of the ten deciles of the top 250 stocks, with each decile comprising 25 stocks. The forward and 
backward flows need not be equal, and they do not appear to be equal in Figure 4. We sch^ from the 
chart that for the largest 250 stocks the flow is downward. For the smaller stocks, we would expect 
the flow to be upward. 

If we follow the flow of a stock that occupies a given rank at time zero, then the expected rank 
of the stock will change over time according to its flow. Suppose a stock is at rank k at time 0, and 
let us estimate its expected rank at time r e K by 

R-A:(t)= Jim 1 / rs+T{Vsi}i))ds. 

In this case, Rfc(O) = k, and if this rank is among the higher ranks, we would expect the flow to be 

negative, which would mean that for r > we would expect that Rfc(T) < k and Rfe(— r) < k. We 
would like to use the to estimate the gfe, and although Rfe(T) need not equal Rfe(— r), the gk 
generated using either one will provide estimates for the solution of (3.6). Accordingly, we shall use 
the average of the two, with 

Rfc(r)+Rfc(-r) 

2 J' 

where the brackets signify the nearest integer. Values of Rfe(T) for fc = 1, . . . , 250 and r ~ ±4 are 
shown in Figure 5, and the values for positive and negative r are clearly different. We have no 
explanation for this difference. 

Figure 4 was generated by following the market weights of stocks that occupied a given rank at 
a given time in the decade from January 1, 1990 to December 31, 1999. Since stocks enter and leave 
the market, we used only the largest 250 stocks, after eliminating any stocks that did not have a 
full ten-year history. The trajectories of the weights for each to the top 250 ranks were followed for 
1000 days forward or backward, and were then averaged over all starting dates that would allow 
the full 1000 days to be used. Finally, the ranks were separated into deciles, with ranks 1-25 in the 
first decile, 26 50 in the second decile, and so forth. The curves in Figure 4 represent the average 
trajectories for the weights of each of the ten deciles, forward and backward. 

Figure 5 was generated by following the weight trajectories used for Figure 4 and, for each 
trajectory, noting the starting rank and ending rank, i.e., the rank after 1000 days (approximately 
four years of trading days). The final rank corresponding to the initial rank k in Figure 5 is the 
average ending rank for those trajectories that begin at rank k at time 0. This was carried out in 
forward time and reversed time. 
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Figure 4: M(fc) (0)e'^'=(^) (black, solid), Ai(fc) (0)e'^'=('') (red, dotted): 1990-1999. On average, 
a stock that starts at time at rank k with market weight /i(fc)(0) will move to weight 
A*(fc)(0)e'*'=^^'' at time r £ [0, 1000], or to weight Ai(fc)(0)e'^'=(^' in reversed time. 
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Figure 5: Rfe(4) (black, sohd) and Rfc(-4) (red, dotted): 1990-1999. 
On average, a stock that starts a given initial rank will move to the 
corresponding final rank four years later, or earlier, in reversed time. 
The straight line represents final rank = initial rank. 



10 



Let Gfc(r) be the expected growth rate at time r e K of a stock which occupies rank k at time 
0, and we shall estimate Gfe(T) and Gfe(— r) from the slope of the forward and backward flow at 
rank k, so 

Gkir) =DrMr) and Gk{-T) = 0^^^), (4.4) 
for r > 0, with Gfe(O) = Qk- We shall use the average 

Gfe(T)^^(Gfe(T) + Gfe(-T)) 

to estimate the rank-based growth rates gk- In the data we analyzed, the derivatives in (4.8) at 
r = 4 were estimated by measuring the rate of change of the flows 0(t) and ^(r) for the period from 
day 981 to day 1000, and then annualizing this rate. 

Since for a given stock the name-based growth rate is invariant with rank, the same holds for 
the average of the name-based growth rates weighted by occupation rates, 

n 
i=l 

Hence, 

n 

'^k{T)=9^^(r)+Yykili, (4.5) 

i=l 

for T e M, where 

9K,(r) = {^ + 1 - Rfe(T)) .g, + (Rfe(T) - 

and t the largest integer such that ^ < Rfe(T). If we combine (4.5) with (3.3), we find that 

(r) =9k + Gk{T) - gk. (4.6) 

We can first estimate gk, Gfc(r), and Rfe(T), and then use (4.6) to recursively generate the values 
of the rank-based growth rates gk for a subsequence of ranks of the form k, Rfe(T), Rr;j.(^) (''"), • • •, as 
well as interpolated points. 

Once we have estimates for the values of the gk, we can estimate the -ji directly by using 

7i = ^ ( Jim 1 / {d log Hi {t) - g^, (j) dt) + lim 1 / (d log Jli (t) - (j) dt) ) . (4.7) 

Our second-order model for the market X will then be 

dlogXi{t) = (7i + gft{{))dt + cTf^f^i-) dWi{t). 

The various steps in the estimation process are shown in Figures 6, 7, and 8. In Figure 6 we see 
the estimated forward rank Rfe(4) versus the initial rank k, and find that the relation is quite close 
to linear with 

Rfe(4)^4.6 + 1.16fc. 
With this estimate, we can use (4.6) in the form 

9(4.6+i.i6k)=9k + Gfc(4) - gk (4.8) 
to estimate the 9k from the values of gk and Gfe(4). 
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The values of Gfe were estimated from the slopes of the flows used to generate Figure 4 for the 
ten rank-decile groups of 25 stocks each from the largest 250 stocks. Linear approximations for 
all the ranks were generated using a least squares flt. These results appear in Figure 7, and the 
corresponding linear equations are 

Gfc(O) = -4.2 - .034fc and 6^(4) = -4.5 - .027k. 

By using the values derived from these equations in (4.8) we can generate values for for an 
increasing sequence of ranks k. The chart in Figure 8 shows these values with linear interpolation 
connecting the points to generate a continuous curve. 
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Figure 6: Estimated four-year forward rank Rfe(4) corresponding to initial rank k. 

Rfe(4)^4.6 + 1.16A;. 
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Figure 7: Estimated expected growth rates Gfc(O) and Gk{4:) at rank k. 
Gk{0)^ - 4.2 - .034fc (black, solid line; dots), 
Gfe(4)= - 4.5 - .027fc (red, broken line; circles). 
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Figure 8: Values of gk for ranks 1 to 250, calculated recursively and interpolated from 
9(4.6+1. 16k)— dk + Gfc(4) — Qk- The gk here are not normalized to add up to 0. 
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Once we have estimates for the values of the gk, we can use these values along with (4.7) to 
estimate values for the 7^ of individual stocks by name. The integrals in (4.7) were approximated by 
daily logarithmic relative returns taken for each of the stocks along with the values of the Qk- The 
(non-normalized) values for 7j for the decade 1990-1999 for a number of well-known stocks appear in 
Table 1. This is hardly a definitive study, so only a few stocks are included here. Moreover, in some 
future work, it would be desirable to have confidence intervals for these values, rather than point 
estimates. In that regard, probably the most promising method would use some form of jackknife 
estimator, with perhaps 12 pseudovalues generated by leaving out one month of the year at a time 
(see Mosteller and Tukey (1977)). Probably the entire estimation process would need to be repeated 
for each pseudovalue. 

The values in Table 1 were estimated using combined forward and backward estimates, as in 

(4.7), for the decade 1990-1999. Using only forward estimates or only backward estimates for the 
7j could have produced biased estimates, since some of these companies grew considerably over that 
decade. For each company, the number in parentheses is the rank of the time-averaged log-weight 
of the stock during the decade. 

While the values in Table 1 may not be definitive, they at least appear plausible. The higher- 
ranked stocks have generally higher 7, which should help them maintain their positions at the top 
of the market. At this writing, Apple, AAPL, has the highest market capitalization in the U.S. 
market, but in the 1990s we see that its average rank was 93, and its 7 is correspondingly —1.67%. 
Hence, the estimated 7^ provide no miraculous forecasts of future behavior; instead they reflect local 
stability consistent with the observed decade. 



5 Conclusion 

The purpose of first- and second-order models for stock markets is to create a rigorous backdrop 
for the statistical analysis of the behavior of individual stocks. Second-order models provide a 
more accurate and complete representation of a stock market than is possible in first-order models. 
The estimation of parameters for second-order models is more involved than for first-order models, 
and implicit methods must be used. We have proposed methods for the estimation of second- 
order growth rate parameters, and with these methods a more complete stock-market model is 
possible. Nevertheless, our techniques are rudimentary, and we believe that future research will 
yield significant improvements. 



Table 1: Values of 7^ for various companies, 1990-1999. 



Apple, AAPL (93) 
Coca Cola, KO (4) 
Exxon, XON (3) 



1.67% 
0.26% 
0.11% 
0.14% 
0.10% 
0.12% 



General Electric, GE (1) 

International Business Machines, IBM (6) 

Microsoft, MSFT (5) 
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